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SUMMARY 

In order to predict the performance of a synthetic aperture radar (SAR) for 
a wide variety of system and environmental parameters a computer simulation 
of the SAR system has been developed. This model is a detailed description 
of the SAR imaging process on a pulse-by-pulse basis. The simulation is im- 
plemented as seven computer programs fcr a CDC Cyber 171 digital computer. 

The first of these programs is the simulation controller. This program co- 
ordinates communication between the other programs and between the user and 
the simulation. It is also responsible for computing certain parameters, 
like the map start time, that will be needed by the other programs. The 
second program discussed in this paper specifies the model of the terrain to 
be mapped. The terrain can consist of discrete scatterers or extended homo- 
geneous areas. The homogeneous areas are represented by a collection of 
many discrete targets. The density of the discrete targets per unit area 
is under user control but the scatterers are placed so as to insure several 
per SAR resolution cell. Representing homogeneous areas as collections of 
point reflectors enables the simulation to support the phenomenon of co- 
herent speckle which is characteristic of some monochromatic SAR imagery. 

The third program discussed computes the range between each scatterer in the 
terrain and the phase center of the antenna. These calculations are made on 
a pulse-by-pulse basis for each pulse in the simulated time interval. The 
range is a function of the location of the terrain on the planet surface, 
the rotation and shape of the planet, the orbital parameters of the satellite 
carrying the SAR and the location of the antenna phase center with respect 
to the center of mass of the satellite. This program also weights the re- 
ceived echo strength in accordance with the antenna gain pattern. The 
fourth program computes the IPQ (In Phase arid Quadrature) video signal by 
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convolving the radar impulse-response and the echo time histories. The 
fifth program handles motion compensation. It applies the Doppler phase 
corrections for each pulse to effect SAR beam steering and range focusing. 

The sixth program is the SAR processor. It weights the IPQ video and com- 
putes the fourier transform to produce the complex output image. The seventh 
and last program is the post processor# It computes the image log magnitude 
leveJs, assigns grey shade levels and displays the resultant image. The 
output of the simulation is a fully processed SAR image displayed on a high 
resolution CRT. 

1.0 PURPOSE AND STATUS 

In general the simulation is designed to permit investigators to assess the 
effect of various radar and geometric parameters on the quality of the SAR 
imagery produced. In particular it will be used to investigate the follow- 
ing: 

1. To assess the effect of various antenna patterns on the SAR image, 

2. To assess the effect of various SAR parameters on the system's 
ability to detect small differences in the radar cross section of 
adjacent homogeneous areas. 

3. To assess the effect of mapping geometry on SAR image quality. 

The first SAR simulation produced at ARL:UT modeled an airborne SAR scenario. 
This simulation has been completed and results obtained from this first 
simulation indicate that our approach was feasible. The simulation of the 
orbital case, on a faster computer, is nearing completion. 

2.0 PROGRAM STRUCTURE 

The structure of the simulation is shown in Figure 1. The simulation con- 
trol program, SARTREK, is run first. It reads all the user input data, com- 
putes various system parameters based on these inputs and stores these data 
in an area which is accessible to the other programs. As the other programs 
in the simulation are executed they will produce other parameters necessary 
to subsequent programs. These data are also stored in the "control para- 
meter" area. The user supplies input data which specifies the antenna 
pattern, terrain parameters, SAR radar and processor parameters, orbital 
parameters and a description of the planet being mapped. After SARTREK the 
data flow is from left to right in Figure 1. Each program is run in turn 
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FIGURE I 

PROGRAM STRUCTURE FOR SAR SIMULATION 














and its output is stored on tape or disk. The output of one program forms 
the inpuft to the. next. This technique has two advantages. First, ic re- 
duces the amount, of computer time one must use in a single block. Second, 
it yields a considerable savings in computer time since the investigator can 
change parameters in a given program and rerun the simulation without having 
to rerun the program upstream from where the chang( « were made. For example, 
the user can investigate the effect of changes in Doppler processor para- 
meters (SARP) without rerunning the four programs which precede it. The 
output of the program is a fully processed SAR map of the modeled terrain 
displayed on a high resolution CRT. 

3.0 SARTREK - THE SIMULATION CONTROLLER 

Program SARTREK forms the communication link between the user and the simula- 
tion and between the various programs within the simulation itself. SARTREK 
performs four major functions; it reads the user specifications of the 
simulation to be run, it computes the mapping start time for the conditions 
specified; it computes the change in incident angle during the formation of 
the SAR map and it computes the range of useable pulse repetition frequen- 
cies (prf). Each of these functions will be discussed in turn. 

Program SARTREK reads the specification of the simulation parameters su - 
plied by the user. These spt cifications are listed below. The user speci- 
fies certain parameters describing the planet about which the SAR is 
orbiting. These parameters are; 


1. 

Planet 

Name, 

2. 

Planet 

Equatorial Radius, 

3. 

Planet 

Eccentricity, 

k. 

Planet 

Gravitational Constant , 

5. 

Planet 

Rotational Rate and 

6. 

Time of Prime Meridian Passage. 


The planet name is used simply as an identifier for the set of data which 
describes the planet. The pquatorial radii , and the planet eccentricity 
determine the shape of the planet. The length of the radius vector from 
the center of the planet to a po:nt on it’s surface is a function of lati- 
tude. The planet is assumed to be symmetric in longitude. By modeling the 
planet, in this manner one can investigate the degradation in a SAR map due 
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to assuming it to be a spherical planet. The planet gravitational constant 
is necessary to define the satellite orbit. The planet rotational rate and 
radius at the point being mapped are necessary ingredients in the computa- 
tion of the distance between scatterers in the terrain and the antenna phase 
center on a pulse-to-pulse basis. Unlike the aircraft case, uhe terrain 
cannot be considered stationary in the orbital situation. Planet motion 
influences the Doppler shift observed from the terrain and introduces range 
gate positioning problems. The time of prime meridian passage is the time 
at which the planer's prime me rid.* an crosses a fixed reference point in 
inertial space. For the earth this reference is the point at which the 
vernal equinox occurs. 

The user specifies the orbital parameters to be us^*d. They are listed below 

1* Orbit Identification, 

2. The Semi-Major axis of the orbit, 

3. The Orbit Eccencricity, 

1+. The Inclination of the Orbit, 

5. The Argument of the Perigee, 

6. The Longitude of the Ascending Node and 

7. The Time of Perigee Passage. 

Again the orbit identification is a serial number used to identify the set 
of parameters with which it is associated. The semi-major axis of the orbit 
ellipse determines the scale of the orbit. The orbit eccentricity deter- 
mines the shape of the orbit. An eccentricity of 0 specifies a circular 
orbit and an eccentricity of one specifies a parabola. Eccentricities bet- 
ween zero and one yield an elleptical orbit. The inclination of the orbit 
specifies the angle between the plane of the orbit and the equatorial plane 
of the planet being orbited. The longitude of the ascending node specifies 
the orientation of the orbit i\l plane about the north-south axis of the 
planet. The argument of the perigee specifies the orientation of the major 
axis of the orbit ellip ^ in the plane of the orbit. Lastly, the time of 
perigee passage determines the location of the radar platform in the speci- 
fied orbit at a given time. 

The user specifies uht parameters of the terrain model. These include; 
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1. The latitude of the terrain, 

2. The number, cross section ana location of discrete scatterers, 

3. The number, size and location of each homogeneous fields and 

h. The scatterer density and strength within each homogeneous field. 
These parameters will be discussed in detail in the section which desciibts 
the terrain model. 

The user also cpeciiMes the SAB radar and processor parameters. The para- 
meters are listed below; 

1. Radar configuration identification, 

2. Radar wavelength, 

3. Sroand range sample interval, 

k. Azimuth resolution, 

5. Azimuth sampling interval, 

6. Radar impulse response shape, 

7. Doppler weighting functions, 

8. Pulse repetition frequency and 

9. Overlay ratio. 

The radar configurat .on identifier is a serial number used to denote the set 
of parameters with which it is associate' 1 Most of these parameters need no 
additional exnlanation. The user specifiers the ground range sample interval 
desired. Once SARTREK has determined the mapping geometry it determines 
and reports to the user the slant range sample interval r. quired to yield 
the desired results. The radar impulse response shape, which determines the 
system range resolution, i s specif, ed in terms of the slant ran - e sample 
interval. 

The user specifies the mapping geometry by indicating the desired nadir and 
squint angle and the swath width. 

Tne radar platform is specified by the orientation of the platform body axis 
at the start of mapping and the platform roll, pitch and yaw rates. 

The antenna is specified by the following; 

l. Antenna identification, 

2. Antenna gain pattern, 
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• Location of the ante^ u phase center with respect to the 
platform body axis coordinate system* and 
U. Orientation of the ante.ina coordinate system with 
respect to the platform body axis coordinate system. 

Again* the antenna identification is just a number used to identify a parti- 
cular set of antenna parameters. Hie antenna gain pattern is read as a 
table of values which specify antenna gain as a function of the bores ight 
angle through which the energy passes to and from the antenna. By permit- 
ting the user ;o separate the phase center of the antenna and the center of 
mass of the platform* the simulatior can be used to model the effects of im- 
proper stabilization of the radar platform on the image. Similarly, the 
orientation of the antenna coordinate system can be used to study the ef- 
fects of antenna pointing errors. 

Program SARTREK also c es the map start time and the lergth of the SAR 
array. Hie user has spr.ifted the platform orbit, the planet, the latitude 
at which the terrain to be mapped is located and the mapping geometry. The 
programs takes «nese parameters and computes those times in the satellite 
orbit for which tl specified conditions are met. SARTREK produces the plot 
shown in Figure 2. This plot is used to provide the user with a general 
view of the results of his specifications. The ground track of two orbits 
are shown. In the figure shown the orbit is circular and inclined at 45° 
to the equator. The user has set the terrain at ^5° north latitude, tne 
nadir angle at 60° and the squint angle at 90 ° to the xeft of the velocity 
vector. The program begins plotting the -round track at start of the first 
mapping time. For the case illustrated the ground track starts just over 
noi*chv< stem Turkey. The program also draws a line from the ground track 
to the center of the map. For this case the line is perpendicular to the 
satellite ground track reflecting the 90° squint angle specified. Note 
that the program locates three other mapping locations which satisfy the 
requirements during the two orbits. The simulation uses the first one 
found. 

Program SARTREK also computes and reports to the ;.ser the change in incident 
angle which occune- '.iring the time of map formation. The program computes 
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the length 01 the array necessary to provide the desired resolution and 
froa this and the geometry of the given mapping situation computes the 
change of incident angle during the map. This information is passed to the 
terrain simulation to be used there to change the reflectivity of the ter- 
rain if the change in incident angle is large enough to warrant such a 
change. 

SARTREK also computes a suggested range of prf's to the user. These prf 
ranges are based on the following criteria. First, the prf must be high 
enough to permit at least two pulses per resolution cell. Second, the prf 
should be low enough to preclude range ambiguities within the range extent 
of the antenna main lobe footprint on the gromd. Thirdly, the prf should 
be high enough to preclude azimuth ambiguities within the cross range extent 
of the antenna mainlobe footprint. Fourthly, the prf must preclude eclip- 
sing, that is, the system must not emit a pulse during the time echoes are 
being received from the desired swath. Prf ranges which satisfy these 
requirements are reported to the user. However, the user is free to select 
whatever prf he desires. 

U.O TERRAIN - THE TERRAIN MODEL 

Program TERRAIN produces a set of data which represent the terrain. In 
general terrain can be divided into two types of scatterers, discrete and 
extended. Discrete or point scatterers were usually relatively small ob- 
jects , often man-made such as automobiles , buildings, etc., any object that 
is characterized by strong, specular, localized reflection. The extended 
or homogeneous scatterer is a larger target such as a forest or wheat field. 
The cross section of an extended target is fairly uniform within a given 
area. For many purposes the cross section of an extended target is speci- 
fied as cross section per unit area. This leads to simulating such a target 
as: a bounded area, the interior of which is -»et to some cross section per 
unit area. This representation is not sufficient for representing homo- 
geneous targets in t.ie SAR application. As a result of coherent processing, 
homogeneous areas of a SAR image can exhibit an effect known as coherent 
speckle. Cohere. eckle appears in the SAR image as dark winding trails 
of various lengths. The origin of coherent speckle lies :*n the fact that 
backscatter from a "homo^c neous ,f area is really due to reflection from many 
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scatter points vithin each resolution cell. For a coherent processor the 
received power from each cell depends on the complex sun of the echoes from 
the many discrete scatterers in each resolution cell. Thus two resolution 
cells containing the same umber of identical scatterers would yield dif- 
ferent received power due to differences in the spacing of the scatterers 
vithin the resolution cell. The degree to which coherent speckle affects 
the image is a function of the width of the mainlobe response with respect 
to spacing of the scatterers within the resolution cell, the sampling inter- 
val and the amount and type of incoherent integration. The presence of co- 
herent speckle increases the variance in the measurement of cross section 
in a homogeneous field and thus increases the error in the estimate of cross 
section. Since one purpose of the simulation is to assess the effect of 
various SAR parameters on the system* s ability to detect small differences 
in cross section of adjacent homogeneous areas, it is necessary that the 
simulation support the origin and propagation of coherei t speckle. For this 
reason the terrain program models homogeneous areas as a collection of many 
point scatterers. The user specifies each homogeneous field as a rectan- 
gular area. The center and extent of the area is specified with respect to 
the center of the terrain. Hie user also specifies the strength of the 
scatterers and the number per resolution cell. Program TERRAIN then com- 
putes the required number of scatter ers and places them at random within the 
extent of the field. Hie user may specify up to nine such homogeneous 
areas* Up to fifty discretes may also be specified. Note + bat the loca- 
tion of each scatter is specified with respect to the center of the terrain 
being modeled. Placement of the terrain on the planet surface is accom- 
plished by the next program in the cuain. This permits the terrain to be 
placed at any desired latitude without rerunning program TERRAIN. At the 
user's option TERRAIN produces a plot of the terrain for inspection. 3uch 
a plot, showing a single homogeneous field and four discrete scatterers is 
shown in Figure 3. 

5.0 SLNTRNGE - THE SLANT RANGE COMPUTATION 

Hie program SLNTRNGE has three functions. The first of these is to compute, 
for each pulse, the slant range between each scatterer in the terrain and 
the phase center of the antenna. Figure U indicates the geometry involved. 
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FIGURE 9 
TERRAIN OUTPUT 

HOMOGENEOUS FIELD AND FOUR DISCRETE SCATTERED* 






The location of each scatterer is known with respect to the terrain co- 
ordinate system. *1116 location of the center of the terrain is known 
relative to the planet center since the user has specified the latitude of 
the terrain ana SARTREK has computed the longitude at which it must be 
placed to be mapped from the given orbit and geometry. The latitude and 
longitude yield the direction of the vector from the planet center to the 
terrain center. The radius of the planet at the latitude indicated gives 
the length of the vector. The orbital parameters yield the vector from the 
planet center to the center of mass of the SAR platform. The vector from 
the platform mass center to the phase center of the antenna is dependent 
cm the location of the antenna with respect to the mass center and the 
orientation of the platform body axis. The vector between the scatterer and 
the antenna phase center can be found once the other four vectors are known. 
This computation involves several coordinate transformations repeated for 
each scatterer and each pulse. Much attention has been given to making the 
code that accomplishes this as efficient as possible. 

A y ther function of SLNTRNGE is to compute the angle between line-of-sight 
t. each scatterer and the antenna boresight. The gain to be associated 
with this angle is looked up in the antenna gain pattern table. This gain 
is then applied to the scatterer strength. This also is a potentially time 
onsuming task. Time has been saved in this area by computing the antenna 
angles subtended by the terrain patch for et- v pulse and restricting the 
gain pattern search within these bounds. 

The final function of SLNTRNGE is to reduce t no echo strength in accordance 
with the space loss over the slant range. 

The output of SLNTRNGE is a set of -^ant ranges and their associated echo 
strengths for each scatterer cn each pulse. Note that for each pulse 
SLNTRNGE computes the range to all the scatterers. These ranges change on a 
pulse-to-pulse basis as the satellite moves. It is this range change with 
time that appears as a Doppler frequency in the SAR processor. Program 
SLNTRNGE will plot the range profiles for a given pilse at the user’s option. 
Figure 5 shows two examples of these plots. The upper plot was produced 
from a terrain model which contained seven discrete targets. The lower 
plot was produced from a terrain model consisting of a single homogeneous 
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field. The reduction is density of the scatterers at the near and far edge 
of the terrain is caused by napping the field from a diagonal aspect. Both 
Of these plots were produced by the airborne SAR simulation. 

6.0 CONVOLVE - THE PRODUCT! OH OF THE IPQ VIDEO. 

The objective of Program CONVOLVE is to compute the interaction between the 
radar signal, the terrain and the radar receiver to produce the IPQ video 
signal to be processed by the SAR processing portion of the simulation. The 
overall simulation can be divided into two parts. The first part starts 
with the terrain and progresses to the IPQ video (Programs TERRAIN , SLNTP-NGE 
and CONVOLVE ) . The second part begins with the IPQ video and moves back to 
a representation of the terrain which is the SAR map (Programs MOCO, SARP 
and SARD). 

The anten: w and radar are assisted to be linear networks. Under this assunp- 
tion the output of the radar receiver is the convolution of the received 
signal and the radar impulse response, 

H(x) = F(x) • G(x), (1) 

where 

F(x) is the input signal in the time domain, 
u(x) is the radar impulse response function, 

H(x) is the IPQ video and 
* is the convolution operator. 

The input signal, F(x), has been computed by SLNTRNGE and the radar impulse 
response function has been specified by the user. The convolution can be 
computed directly or by using Fourier transform theory. Recall that the 
convolution of two functions is the transform of the product of their trans- 
forms. Using bars to denote Fourier transforms. 


H(x) = F(x) • G(x) (2) 

This formulation greatly speeds the convolution process. The convolution 
must he computed once for each pulse. The impulse response function does 
not vary with time so G(x) need only be computed once. Figure 6 shows the 
result of this process. Diagram a shows the impulse response and diagram b 
shows the output of SLNTRNGE. Diagrams c and d are the I and Q components 

i 

*. i * 

*■ .;u 
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of the IPQ video obtained from the convolution. 

7.0 M0C0 - MOTTOS COMPENSATION 

Program MOCO has two objectives. The first is to introduce the effects of 
motion compensation or focusing on the IPQ video waveform. The second ob- 
jective is tc accomplish the orthogonal read or "memory corner turning" 
necessary to further processing of the SAR signal. Focusing the synthetic 
array on a point in space is accomplished by adjusting the phase of the IPQ 
signal for each pulse so that thr total time required for a signal to travel 
from the focus point to the phase center of the antenna is the same for all 
pulses. In the simulation the location of the scatterer with respect to 
the antenna phase center is known exactly. Thus perfect motion compensation 
or focusing is possible. In real systems this is .. * the case. Errors in 
the location of the antenna phase center with respect to the map center can 
be introduced into MOCO to study the effects of errors in motion compensa- 
tion. 

The second function of MOCO is the orthogonal reordering of the data. SAR 
data are accumulated on a pulse-to-pulse basis. Each pulse contains a 
sample at each range interval. The Doppler processor needs the data pre- 
sented in range bin order. That is, the processor needs the time histories 
of each range bin. 

8.0 SARP - THE SAR PROCESSOR 

In one sense of the word, this portion of the SAR model is not a simulation 
at all. The digital Doppler processor in this simulation differ from an 
operational SAR processor mainly in the machine on which it is implemented. 
Program SARP i~ a basic Fourier transform processor. The IPQ video from 
MOCO is weighted by the function of the user choice. These data are then 
transformed range bin range bin and resultant spectrum is passed to the 
next program in the simulation. If the user desires, program SARP will 
produce a perspective plot of the filter magnitudes for each range bin. 
Figure 7 shows such a plot for a terrain model consisting of seven discrete 
targets . 
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9.0 SARD - THE SAR DISPLAY 

The purpose of program SARD is to display the filter magnitude data as the 
final SAR map. The program has three basic functions. The first is to over- 
lay the filter magnitudes if the overlay option is selected by the user. If 
overlay is selected the corresponding filter magnitudes from the maps in- 
volved are averaged and the result becomes the filtei magnitudes displayed. 
The purpose of this noncoherent integration is to reduce coherent speckle 
and the effect of sampling grid placement. 

The second function is to convert the filter magnitudes to grey shade levels. 
The high resolution has 16 grey shades. The grey shade assignments in.olve 
mapping the filter magnitude values into grey shades in accordance with the 
user’s specifications. 

'Hie final function of SARD is to communicate the grey shade data from the 
computer t^ the display memory where it is transferred to the CRT. Figure 8 
shows photographs of various modeled terrains. These images were produced 
by the airborne SAR simulation. The radar was illuminating the terrain from 
a near diagonal aspect. Detail "a" shows the SAR image for a terror* model 
consisting of seven discrete scatterers arranged along the front and back 
sides of a square with one scatters in the center. For these images the 
resolution cell was not square, this led to the deformation of the image 
that can be observed In the figure. Derail r b" is the image from a single, 
square homogeneous field. Detail "c" is ^he image for a homogeneous field 
with large discretes at each corner. 
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